Atmospheric Solar Heating in Minor Absorption Bands 


Ming-Dah Chou 
Laboratory for Atmospheres 
NASA/Goddard Space Flight Center 
Greenbelt, Maryland 
USA 


November, 1998 

Terrestrial, Atmospheric and Oceanic Sciences 


Corresponding author: Dr. Ming-Dah Chou 

Code 913 

Laboratory for Atmospheres 

NASA/Goddard Space Flight Center 

Greenbelt, MD 20771 

Telephone: 301-286-4012 

FAX: 301-286-1759 

e-mail chou@climate.gsfc.nac3.gov 



ABSTRACT 


Solar radiation is the primary source of energy driving atmospheric and oceanic 
circulations. Concerned with the huge computing time required for computing radiative 
transfer in weather and climate models, solar heating in minor absorption bands has often 
been neglected. The individual contributions of these minor bands to the atmospheric 
heating is small, but collectively they are not negligible. The solar heating in minor bands 
includes the absorption due to water vapor in the photosynthetically active radiation (PAR) 
spectral region from 14284 cm' 1 to 25000 cm' 1 , the ozone absorption and Rayleigh 
scattering in the near infrared, as well as the 0 2 and C0 2 absorption in a number of weak 
bands. Detailed high spectral- and angular-resolution calculations show that the total effect 
of these minor absorption is to enhance the atmospheric solar heating by -10%. Depending 
upon the strength of the absorption and the overlapping among gaseous absorption, 
different approaches are applied to parameterize these minor absorption. The 
parameterizations are accurate and require little extra time for computing radiative fluxes. 
They have been efficiently implemented in the various atmospheric models at 
NASA/Goddard Space Flight Center, including cloud ensemble, mesoscale, and climate 
models. 
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1. INTRODUCTION 

The ultimate force driving the Earth’s weather and climate is the solar radiation. 
Satellite measurements of the earth radiation budgets show that the Earth reflects -30% of 
the incoming solar radiation, and the rest is absorbed at the surface and in the atmosphere 
(Barkstrom, 1984, Kandel et al., 1998). Partitioning of the 70% of the solar heating 
between the surface and the atmosphere is very important to oceanic and atmospheric 
circulations; a larger solar heating at the surface will cause a higher surface temperature and 
an enhanced evaporation, leading to a more energetic atmospheric circulation. However, 
this partitioning cannot be directly estimated from currently available satellite 
measurements. It can only be estimated either from model calculations or from surface 
measurements with accuracies not well known. The Atmospheric Radiation Measurements 
(ARM) program of the US Department of Energy has conducted a number of field 
experiments to study the radiative processes in the atmosphere. Both narrow- and broad- 
band flux measurements are compared with theoretical model calculations. The aim is to 
provide information for validating radiation parameterizations used in numerical weather 
and climate models. It is clear that accurate calculations of solar heating is important in 
climate studies using computer models. 

The broadband parameterizations for solar radiation used in weather and climate 
studies are mostly based on high spectral resolution calculations of gaseous absorption and 
detailed angular integration of particle scattering. To enhance computing speed, the 
spectrum is divided into a small number of bands, and only important gaseous absorption 
and particle scattering are included. Minor absorption and scattering due to gases or 
scatterers spread over a wide spectral range. They are almost always neglected in radiation 
models used in weather and climate studies in order to save computing time. Individua y 
their effects are small, but the total effect may not be negligible. As the capability of the 
computing system increases, more detailed physical processes are included and higher 
spatial and temporal resolutions are used in weather and climate models. To be consist t 
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with these improvements, radiative transfer calculations are required to be more accurate 
than the existing models. Because calculations of radiative terms in weather and climate 
models require a large portion of the total computing time, it is essential that improvement 
in radiation parameterization should not have a large impact on computing speed. 

There are a number of solar radiation schemes available for use in global climate 
models (cf. Fouquart et al., 1991). At NASA/Goddard Space Flight Center, a series of 
parameterizations for radiative transfer in the solar spectrum due to various absorbers and 
scatterers has been developed and applied to global and regional climate studies (Chou, 
1986; Chou, 1992; Chou and Lee, 1996). The ozone absorption and Rayleigh scattering 

were included in the ultraviolet (UV) and visible regions with the wavenumber v> 14285 
cm' 1 (0.7 (im) but were neglected in the near infrared. On the other hand, the absorption 
due to water vapor was included in the infrared with v< 14285 cm' 1 but were neglected in 

the photosynthetically active radiation (PAR) region between 14285 cm' 1 and 25000 cm' 1 . 
(0.4 pm and 0.7 pm). Oxygen and C0 2 also absorb solar radiation in the infrared. The 
parameterizations only included the absorption in major 0 2 and CO, bands. Absorption in 
minor 0 2 and C0 2 bands, which scatter over wide spectral ranges, was neglected. 
Individually, these weak absorption and scattering neglected in the parameterizations have a 
little effect on solar heating of the atmosphere. Collectively, their effect is not necessarily 
negligible. This study presents the parameterizations for the minor absorption and 
scattering, which have been efficiently implemented into the mesoscale and global 
atmospheric models used at NASA/Goddard Space Flight Center. 

2. MINOR GASEOUS ABSORPTION AND RAYLEIGH SCATTERING 

In a clear atmosphere, the heating of the atmosphere is primarily due to water vac >r 
and 0 3 and secondarily due to O, and C0 2 . Rayleigh scattering interacts with absorpti n 
and affects solar heating in the atmosphere and at the surface. Both the ozone absorpt' a 
coefficient and the Rayleigh scattering coefficient decrease rapidly with decreasing 
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wavenumber v. Figure 1 shows the spectral distribution of the ozone transmission 
function of the entire atmospheric column in the vertical direction. The Rayleigh scattering 
(not shown in the figure) decreases with decreasing wavenumber according to v 4 . For 

v< 14285 cm' 1 , the ozone absorption and Rayleigh scattering are both weak and were not 

included in the radiation model. Water vapor absorbs solar radiation at all wavelength 
except in the UV region. Figure 2 shows the spectral distribution of the water vapor 
transmission function of the atmosphere in the vertical direction. The width of water vapor 
molecular lines is small compared to the mean line spacing. As a result, the absorption 
coefficient varies rapidly within a narrow spectral interval. The transmission function 
shown in Fig. 2 is the mean over 10 cm' 1 . It can be seen that there are many water vapor 

absorption bands with band-strength decreases rapidly with increasing v. The bands in the 

high frequency region with v>14285 cm' 1 are very weak and were neglected in our 
radiation model. 

In addition to water vapor and 0 3 , oxygen and C0 2 also absorb solar radiation. 
Figures 3 and 4 show, respectively, the spectral distributions of the 0 2 and CO, 
transmission functions of the entire atmospheric column. The two major O, A and B 
bands are located at 13150 cm' 1 and 14550 cm' 1 . There are also minor bands located at 
7890 cm' 1 and 15870 cm' 1 . Only the two major A and B bands were included in our 
radiation model. Carbon dioxide absorbs solar radiation in a wide spectral range. Only 
the absorption in 3600 cm' 1 and 5000 cm' 1 bands were included. The absorption in the 
2325 cm' 1 (4.3 pm) band is very strong, but the insolation (incoming solar radiation) is 
small. Solar heating due to C0 2 in this band is small and was not included in the 
radiation model. There is also absorption due to C0 2 in the spectral region between 6000 
cm' 1 and 7000 cm' 1 (1.43 pm and 1.67 pm). The solar energy contained in this spectr 1 
regions is large, but the absorption coefficient is small. It was also not included in the 



radiation model. Table 1 shows the major absorption and scattering included in our 
radiation model, and Table 2 shows those minor absorption and scattering which were 
not included in the radiation model but are parameterized in this study. 

3. HIGH SPECTRAL-RESOLUTION CALCULATIONS 

Parameterizations for solar fluxes in minor absorption bands are based on high 
spectral- and angular-resolution calculations. The temperature, humidity, and ozone 
profiles are typical of a mid-latitude summer atmosphere taken from Anderson et al., 

( 1986), which has a vertically integrated water vapor amount of 2.93 g cm : and ozone 
amount of 0.3 18 (cm-atm) STP . The C0 2 concentration is set to be 350 ppmv at all heights, 
and the 0 2 mass mixing ratio is fixed at 0.23 15. 

The 0 3 absorption coefficient varies by several orders of magnitude in the UV and 
visible regions but rather smoothly with wavenumbers. Fluxes are computed with a 
spectral resolution of 1 cm' 1 with the 0 3 absorption coefficient interpolated from the spectral 
values given in WMO (1985). The absorption coefficients of water vapor, 0 2 and CO, 
vary rapidly with wavenumber, and the line-by-line method is used to compute the 
absorption coefficients. The 1996-version of the molecular line parameters compiled at 
AFGL (Rothman, 1987) are used in the line-by-line calculations. The shape of an 
absorption line is assumed to follow the Voigt function. Because the shape of the far 
wings of a line is not well known, the absorption coefficient is set to zero at wavenumber 
>10 cm' 1 from the line center. The spectral resolution of the line-by-line calculations is 
chosen to be 0.005 cm' 1 , which is adequate for flux calculations in the troposphere. 
Rayleigh scattering is computed using the discrete-ordinate multiple-scattering algorithm of 
Stamnes, et al. (1988). The number of angular streams used is six. The spectral resolution 
used in computing Rayleigh scattering is either 1 cm' 1 when overlapping with the 0 3 
absorption or 0.005 cm' 1 when overlapping with the water vapor absorption. 

4. BROADBAND TRANSMISSION AND FLUX PARAMETERIZATIC 3 
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4.1 The k -distribution method 

The width of an absorption line of water vapor, 0 2 , and C0 2 is much smaller than 
the average spacing between absorption lines. The width of a line decreases linearly with 
pressure and is only -0.05 cm' 1 in the lower troposphere. As a result, the absorption 
coefficient varies very rapidly within narrow spectral intervals, and accurate radiative 
transfer calculations require a high spectral resolution. For the case of water vapor 
absorption, the high resolution line-by-line method requires flux calculations at >10 6 
spectral points. For an atmospheric layer where temperature and pressure can be assumed 
constant, the wavenumbers with the same absorption coefficient are radiatively identical 
and can be treated as one entity. Fluxes at those wavenumbers with the same absorption 
coefficient need to be computed only once. Therefore, flux calculations can be greatly 
simplified by grouping the wavenumbers with the same absorption coefficient. Within a 
small interval A v where the spectral variation of insolation is small, the integration of 
fluxes over wavenumbers can be replaced by that over the absorption coefficient. The 
mean transmission function of the interval Av can be written as 

r(vv')= j e~ kv ^ dv I Av - ^ e~ kw g(k)dk (1) 
Av 0 i=l 

where vv' is the absorber amount, k v is the absorption coefficient at the wavenumber v, g 
is the k-distribution density function (cf. Arking and Grossman, 1972) such that the 

fraction of spectrum with the absorption coefficient in the between k-^dk and k + ~dk 

is g(k)dk , and a,- (=g(k i )Ak i ) is the k-distribution function. It has been found (Chou and 
Lee, 1996) that the k-distribution method requires only -10 values of k (i.e. n~10) for 
accurate calculations of solar fluxes using the k-distribution method instead of >10 6 sr 'ctrsi 
points using the line-by-line method. 

In the atmosphere where pressure and temperature change with height, the 
wavenumbers with a common absorption coefficient at a given height will not necess; ily 
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have a common absorption coefficient at other heights. These wavenumbers are no longer 
radiatively identical, and the k-distribution method cannot be used without applying 
assumptions. Pressure affects the width of an absorption line. Near the center of a line, 
the absorption is strong, and flux calculations are not sensitive to errors in the absorption 
coefficient. Distant from line centers, the absorption coefficient increases nearly linearly 
with pressure. On the other hand, temperature affects absorption primarily through its 
effect on the line intensity, which is constant for a given line. Thus, the temperature effect 
on absorption is rather smooth with respect to wavenumber. It has been shown by Chou 
(1986) that fluxes can be calculated accurately with the use of the one-parameter scaling of 
the absorption coefficient, 


k v (p,T) = k v (p r ,T r ) i 


( p^ m 

\Pr 


f(T,T r ) 


( 2 ) 


where p r is the reference pressure, T r is the reference temperature, m < 1, and f(T,T r ) is 
the temperature scaling function. As explained in Chou (1986), the optical path between 
the top of the atmosphere and the middle and lower troposphere is long near the center of 
absorption lines. Solar fluxes in middle and lower troposphere are primarily attributable to 
the spectral intervals between absorption lines, and the wing-approximation of (2) can be 
applied to accurately compute fluxes. In the stratosphere, the solar heating due to water 
vapor, 0 2 , and C0 2 is small, and underestimation of the absorption coefficient induced by 
the use of (2) does not have a significant effect on the heating rate. In the upper 
troposphere, the absorption is neither strong nor weak, and the error induced by the scaling 
is expected to have a large impact on flux calculations. However, if we choose p r to be 
the upper troposphere pressure, then p I p r ~ 1 in the upper troposphere, and the scaling 
introduces only a small error in the absorption coefficient. It has been found (Chou, 1936) 
that the solar fluxes can be computed accurately by choosing p r = 300 hPa and m- 0.8. 

Flux calculations are not sensitivity to m for 0.7< m< 0.9. The effect of temperature on : 


I 
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absorption of solar radiation is weak. For the absorption due to water vapor, the 
temperature effect on absorption is approximated by using 

f(T,T r ) = l + 0.00135(7- T r ) (3) 

where 7' r =240 K. For the absorption due to 0 2 , and C0 2 , the temperature effect is 
neglected, i.e. f{T,T r )=l. 

Using the scaling of (2), Equation (1) becomes 


T{w) = Ya l e~ k ‘ iPr ’ Tr)W 

i=i 


where 


(4) 


w(p, T ) = w' 



\Pr ) 


f(T,T r ) 


(5) 


Thus, the scaling of the absorption coefficient is reduced to the scaling of the absorber 
amount, and flux calculations are greatly simplified. It is noted that the basis for the k- 
distribution method with the one-parameter pressure and temperature scaling is the same as 
the correlated k-distribution method (Goody et al., 1989; Lacis and Oinas, 1991; Fu and 
Liou, 1992), i.e. the dominant effect of line wings on the absorption. The accuracies of 
these two approaches are comparable, but the former is much simpler than the latter. 

4.2 Flux calculations 

For a spectral band with a width Av, the direct solar flux at a given pressure level p 
can be written as 


F(w)= je~ kvW S v dv 

Av 


( 6 ) 


where S v is the extraterrestrial solar flux density at the wavenumber v multiplied by the 


cosine of the solar zenith angle, k is the absorption coefficient, and w is the scaled 
absorber amount above the pressure level p in the direction of the solar beam. The mean 
transmission of a band is given by 
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r(w)= je k ^S v dv/ Js v rfv 


(7) 


Av Av 

If we divide a spectral band into m small intervals and apply the k -distribution 
method, Equation (6) becomes 


—k, vv 


m 


n 

m 

F(w)=J j S j 

\e~ k ' w dv 

II 

M 

2><\7 5 7 Av 7 

7=1 

Avj 

i = 1 

-7 =1 


(8) 


where Sj is the mean flux density in the interval AVj, a t j is the it -distribution function in 
the interval A Vj, n is the number of the absorption coefficient k, and 


m 

Av=J j Av j 

7=1 


(9) 


By defining 


s=Y. s i Av i 

7=1 


( 10 ) 


and 


m 


hi = X a i,j S j Av j 1 S 
7=1 

Equation (8) becomes 
F(w ) = ST(w) 

where t(vv) is the mean transmission function of the band Av given by 


(ID 


( 12 ) 


t(w) = /i,<? 


—kj\v 


(13) 


i=l 


and is the flux-weighted k-distribution function at p r -~ 300 hPa and T r = 240 K whir 1 .; 
has the property 

2 > = 1 14 ) 

i=l 
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To simplify flux calculations involving scattering, a single value of the scattering 
coefficient is used for each band. Because scattering of light varies slowly with 
wavenumber, this approximation will not introduce large errors in flux calculations. Using 

the & -distribution method for absorption and the 5-Eddington approximation (Joseph et al., 

1976) for scattering, transmission and reflection functions of each layer are first computed 
for each value of k. The two -stream adding method (cf. Chou, 1992) is then used to 
compute the fluxes components F t - associated with k it i=l, 2, ... , n. The total flux of a 
band is the sum of F, weighted by the k-distribution function, 

n 

F(w) = ^ j h i F i (w) (15) 

t=l 

This approach is applied to computing fluxes due to water vapor absorption in the infrared. 

For the case that either the range of k or the value of k itself is small within a spectral 

band, the transmission can be approximated by $ 

1 

t ( W )= e ~~ kw (16) 

I 

where k is the mean absorption coefficient of a spectral band. This approach is applied to 
computing fluxes due to 0 3 absorption. 

Absorption due to 0 2 and C0 2 are treated differently from the absorption due to water 
vapor and 0 3 . The major absorption due to 0 2 is not weak but occurs in narrow spectral 
intervals. An effective absorption coefficient of 0 2 is derived by fitting the mean 
transmission function r(w) computed from the line-by line method, 

r(w) = e“* Vw (17) 

The 0 2 absorption bands are located between water vapor absorption bands. Thus, water 
vapor has little effect on the solar heating due to 0 2 . Instead of computing fluxes in the 0 2 
bands, the reduction in fluxes due to O, absorption is computed from 
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AF(w) = 5(l-e-^) (18) 

Here S is the insolation at the top of the atmosphere in the 0 2 bands. Because scattering 
algorithms are developed for monochromatic radiation, the band-averaged transmittance 
( 17) can only be applied to clear-skies but not to cloudy skies. 

The absorption due to C0 2 is small but overlaps significantly with the absorption due 
to water vapor. Similar to (6) and (12), fluxes in the C0 2 bands are given 

F(w,u ) = J e~^ kvW+ ^ vU ^ S v dv = Sr(w,u) (19) 

where w and u are the scaled water vapor and C0 2 amounts, k and y are the absorption 
coefficients of water vapor and C0 2 , x (w,u) is the transmittance defined by 

n m . . 

/ x v X"" 1 7 —{k.W+YjU) 

T(w,m) = 2, 2A/ rj (20) 

<=iy=i 

and j is the flux- weighted k -distribution function with the absorption coefficients k t and 
Yj for the two absorbers, respectively. Equation (20) requires (mx n) sets of flux 

calculations, which becomes computationally very expensive when used in weather and 
climate models. Other approaches than the k -distribution method have to be used for 
efficient flux calculations when overlapping of absorption is involved. 

The flux-weighted transmittance involving two absorbers can be pre-computed as a 
function of w and u , and fluxes can be computed from ( 19), or alternatively from 

F(w,u) = F(w)~ AF(w,u) (21) 

where A F(w,w) is the flux reduction due to C0 2 with amount u, 

A F(w,u) = J e~ k ’ w (l - e~ r " u )S v dv (22) 

and the integration is over the entire solar spectrum. A two-dimensional table for A F K w,u) 
is pre-computed using the line-by-line method. 
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Table 1 shows the spectral bands of the solar radiation model used at NASA/Goddard 
Space Flight Center. There are eight bands in the ultraviolet and visible region (v> 14285 

cm' 1 ) and three bands in the infrared region (v<14285 cm '). Also shown in the table are 

the absorbers and scatterers included in the calculation of solar fluxes in each band, as well 
as the equations used for computing the transmission function. The first eight bands 
involves the ozone absorption and Rayleigh scattering. These bands are narrow, and an 
effective ozone absorption coefficient and an effective Rayleigh scattering coefficient are 
used for each band. The effective Rayleigh scattering coefficient is derived by averaging 
the flux-weighted scattering coefficient over a band. For the ozone absorption, we first 
compute r(w) from (7) with a high spectral resolution. The effective absorption 

coefficient, k, is then derived by fitting r(w) with (16). The absorption due to water 
vapor is included in the three infrared bands (Bands 9, 10, and 1 1). The water vapor 
absorption varies strongly within small spectral intervals, and the k -distribution function 
method, Equations (13) and (15), is used to compute fluxes. The number of the absorption 
coefficient used in each band is 10. The use of this number is to avoid oscillation in the 
vertical profile of the computed heating rate when a smaller number is used (Chou and Lee, 
1996). The absorption due to 0 2 and CO, occur in wide spectral ranges. The flux 
reduction due to 0 2 is computed from (17) and that due to CO, is computed from A F(w,u) 
using table look-up. Because clouds have a large impact on radiation, we can neglect the 
weak absorption due to O, and CO, within clouds. The flux reduction A F is computed 
only for clear skies or the level above the cloud top. Below the cloud top, it is set to that at 
the cloud top. 

5. ABSORPTION IN MINOR BANDS 

5.1 Water vapor absorption in v > 14285 cm' 1 

The absorption due to water vapor in the photosynthetically active radiation (P/ .} 
region between 14285 cm' 1 and 22700 cm' 1 (0.7 pm and 0.44 |im) was not included ; 



12 


Band 8 of our radiation model. Because the absorption due to water vapor is weak in this 
spectral region, it can be represented by a single effective absorption coefficient k in (16). 
To derive k, the mean transmission function r(vv) given by (7) is first computed using the 
line-by-line method for w ranging from 5 g cm 2 to 20 g cm' 2 . For each w, an effective 
absorption coefficient is computed from 

k(w) = — -lnT(w) (23) 

vv 

The value of k at p r = 300 hPa and T r =240 K is found to range from 0.00065 to 0.00080 
g' 1 cm 2 , and a mean value of 0.00075 g 1 cm 2 is adopted in the broadband flux calculations. 
For the midlatitude summer atmosphere, the vertical-integrated scaled water vapor amount 
in the direction 60° from the zenith is 14 g cm' 2 , and the absorptance in the PAR spectral 
region is -0.0 1 . With a fraction of 0.391 of the extraterrestrial solar radiation contained in 
the PAR, it corresponds to an atmosphere solar heating of -2.6 W m' 2 . Table 3 shows that 
the surface flux reduction in Band 8 due to water vapor is 2.53 W m" 2 from line-by-line 
calculations and 2.42 W m' 2 from the parameterization. Rayleigh scattering and 0 3 
absorption have little effect on the absorption due to water vapor in the PAR. 

5.2 Ozone absorption in v < 14285 cm' 1 

As shown in Figure 1, the absorption due to ozone in the near infrared (Band 9) is 
restricted to a narrow spectral region next to Band 8 where absorption due to water vapor is 
weak. However, the absorption due to water vapor in the rest of Band 9 is not necessarily 
weak, and calculations of ozone heating in the near infrared would require adding another 
band in the broadband radiation model. Because the absorption due to water vapor and 0 3 
is weak in Band 8, the 0 3 absorption in Band 9 can be folded into the absorption in B~nd 8 
as if it were the absorption due to another absorber. To do so, the effective 0 3 absor;. on 
coefficient of Band 8 is enhanced by A k which satisfies 
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(l-<T Mlv ') \s v dv= \[\-e k '' w y v dv (24) 

AV( Av 2 

so that the absorption of solar radiation due to 0 3 in the near infrared is correctly computed, 
where w’ is the ozone amount, and AVj and kVj are the widths of Bands 8 and 9, 
respectively. It is noted that the absorption due to ozone is nearly independent of pressure 
and temperature, it is not necessary to scale the ozone amount given by (5). For a wide 
range of w 1 found in the atmosphere, the value of A k ranges between 0.0032 and 0.0033 
(cm-atm)^' 1 . Therefore, the ozone absorption coefficient in Band 8 is enhanced by 
0.0033 (cm-atm) STP 1 to take into account the absorption by ozone in the infrared. This 
approach to computing the ozone absorption in the near infrared requires no extra 
computing time. For the midlatitude summer atmosphere, the column ozone amount is 
0.3 18 (cm-atm) STP . It can be easily shown that the solar heating of the atmospheric column 
is enhanced by 0.56 W m' 2 , which is independent of the solar zenith angle due to the 
opposite impacts of the solar zenith angle on the insolation and the ozone pathlength. Table 
3 shows that the surface flux reduction in the near infrared due to 0 3 is 0.54 W m' 2 as 
computed from the parameterization, which is very close to that computed with a high 
spectral resolution. 

5.3 Oxygen absorption 

As shown in Figure 3, the absorption due to 0 2 occurs in narrow spectral intervals, 
but is not necessary weak near band centers. In the parameterization of the 0 2 absorption, 
Chou (1990) computed the mean transmission function in the O, A and B bands centered at 
13150 cm' 1 and 14510 cm' 1 and fit the transmission function by (17). The absorption in the 
weak bands centered at 7890 cm' 1 and 15870 cm' 1 was not included. Line-by-line 
calculations show that the total atmospheric heating due to 0 2 is 4.29 W m' 2 for a solar 
zenith angle of 60°. The 0 2 A and B bands contribute 3.70 W m' 2 to the total absorption. 

A heating of 0.49 W m' 2 is attributable to the band 7700-8050 cm' 1 and 0. 1 W m 2 is 
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attributable to the band 15750-15950 cm' 1 . Thus, excluding the two weak bands will cause 
an underestimation of the atmospheric heating by 0.50 W m' 2 . 

To include the absorption in all those bands, the mean transmission function of 0 2 at 
p r =300 hPa and T r = 240 K in the spectral regions 7600-8050, 12850-13190, 14310- 
14590, and 15730-15930 cm' 1 is computed from (7) using the line-by-line method, and the 
effective mean absorption coefficient k is computed from 


k(w) = - 



ln T(vv) 


(25) 


The value of k is between 0.000135 (cm-atm) STP ,/2 and 0.000155 (cm-atm) STP l/2 for a 
large range of w encountered in the atmosphere, and a mean value of 0.000145 (cm- 
atm) STP 1/2 is adopted to compute the flux reduction due to O z from (18). The insolation, S, 
in the spectral regions 7600-8050, 12850-13190, 143 10-14590, and 15730-15930 cm' 1 is 
86.53 W m' 2 or 6.33% of the total solar flux at the top of the atmosphere. Table 3 shows 
that, for a solar zenith angle of 60°, the surface flux reduction due to 0 2 is 4.29 W m 2 from 
line-by-line calculations and 4.18 W m 2 from the parameterization. 

5.4 Carbon dioxide absorption 

The absorption of solar radiation due to C0 2 spreads over the middle infrared from 
2000 to 7000 cm' 1 . It overlaps substantially with the absorption due to water vapor in 
some spectral regions. Line-by line calculations show that, for a solar zenith angle of 60° 
and a C0 2 concentration of 350 ppmv, the atmosphere heating due to C0 2 alone is 8.73 W 
m' 2 but reduces to 3.30 W m 2 when overlaps with water vapor absorption in the mid- 
latitude summer atmosphere. Because of this strong overlapping, the C0 2 transmission 
function cannot be computed independently of water vapor. 

From the pre -computed table of A F(w,u) given by (22), Chou (1990) compute 
flux reduction due to CO, only for the spectral bands centered at -3700 cm' 1 and 5000 on’. 
To include the C0 2 absorption in the minor bands, the flux reduction A F(w,u) is re- 
computed in this study to cover the entire solar spectrum which includes the strong b 
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2480 cm' 1 and the weak bands at 6300 cm' 1 and 7000 cm' 1 . Line-by-line calculations show 
that, for the mid-latitude summer atmosphere, a solar zenith angle of 60°, and a C0 2 
concentration of 350 ppmv, the surface flux reduction due to C0 2 increases by 1.21 W m' 2 
when the minor absorption is included. Table 3 shows that the flux reduction computed 
from the parameterization is 3.60 W m' 2 which is slightly larger than that from line-by-line 
calculations. 

5.5 Ravleigh scattering in the near infrared 

Rayleigh scattering in the near infrared is weak and was neglected in our radiation 
model. To include the Rayleigh scattering in the near infrared, a mean extinction coefficient 
for each of two near infrared bands, Bands 9 and 10, is computed from 

ct = J a v S v dv / j S v dv ( 26 ) 

where the integration is over the spectral interval of a band. The value of o is 0.0000156 
hPa' 1 for Band 9 and 0.0000017 hPa' 1 for Band 10. Rayleigh scattering in Band 1 1 is 
negligible. The effect of Rayleigh scattering in the midlatitude summer atmosphere with a 
solar zenith angle of 60° and a surface albedo of zero is shown in Table 3. Rayleigh 
scattering reduces the surface radiation by 3.12 W m' 2 from the line-by-line and discrete- 
ordinate calculations and by 2.82 W m' 2 using a mean extinction coefficient. 

6. CONCLUSIONS 

Calculations of radiative heating/cooling in the atmosphere require a large portion of 
the total computing time in global climate simulations using a general circulation model. To 
enhance the efficiency of radiation calculations, much of the atmospheric solar heating in 
minor bands are neglected. Individually, these minor heating terms are small, but 
collectively they are not negligible. It is found from high spectral- and angular-resoluti n 
calculations that the absorption due to water vapor in the photosynthetically active radiation 
(PAR) band and that due to C0 2 in middle infrared amount to ~5.5 W m' 2 for a midlat ,;de 
summer atmosphere and a solar zenith angle of 60°. The absorption due to O, and 0 3 : the 
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near infrared amount to -5.0 W m' 2 . Raleigh scattering in the near infrared reduces the 
surface solar heating by -3 W m' 2 . The sum of these heating is -10% of the total solar 
heating of the Earth’s surface, which is not negligible. In this study, different 
parameterizations are applied to these minor absorption and scattering depending upon the 
strength of the extinction and the overlapping of absorption among various absorbers. The 
parameterizations are accurate and require little extra computing time. Therefore, they are 
suitable for climate studies using an atmospheric general circulation model. 
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Table 1. Gaseous absorption and Rayleigh scattering in the broadband 

radiation parameterizations and the equations used for transmittance 
and flux calculations. 


Band 

Spectral Range 
(cm' 1 ) 

Absorber/ 

Scatterer 

Equation 

1 

(44440-57140) 

o 3 

16 



Rayleigh 

-- 

2 

(40820-44440 

o 3 

16 


35700-38460) 

Rayleigh 

— 

3 

(38460-40820) 

o 3 

16 



Rayleigh 

-- 

4 

(33900-35700) 

o 3 

16 



Rayleigh 

— 

5 

(33330-33900) 

o 3 

16 



Rayleigh 

-- 

6 

(31250-33330) 

o 3 

16 



Rayleigh 


7 

(25000-31250) 

o 3 

16 



Rayleigh 

— 

8 

(14280-25000) 

o 3 

16 



Rayleigh 

— 

9 

(8200-14280) 

H 2 0 

13, 15 

10 

(4400-8200) 

h 2 o 

13, 15 

11 

(1000-4400) 

h 2 o 

13, 15 

Total Spectrum 

0 2 

18 

Total Spectrum 

co 2 

21 
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Table 2. Absorption and scattering in minor bands parameterized in this study. 


Band 

Spectral Range 

(cm' 1 ) 

Absorber/ 

Scatterer 

Equation 

Remark 

8 

(14280-25000) 

H 2 0 

16 


9 

(8200-14280) 

o 3 

Rayleigh 

16 

0 3 absorption folded into Band 8 

10 

(4400-8200) 

Rayleigh 

— 


Total Spectrum 

0 2 

18 

New parameterization to 
include minor absorption bands 

Total Spectrum 

CN| 

o 

U 

21 

New parameterization to 
include minor absorption bands 


Table 3. Effects of minor absorption and scattering on the solar heating of the surface 
from detailed high spectral-resolution calculations and parameterizations. 
Fluxes are computed for a typical mid-latitude summer atmosphere. The solar 
zenith angle and the surface albedo are set to 60° and zero, respectively. The 
units of heating are W m' 2 . 


Detailed Calculations 

Parameterization 

Water vapor in PAR 

-2.53 

-2.42 

Ozone in near infrared 

-0.56 

-0.54 

Oxygen in near infrared 

-4.29 

-4.18 

C0 2 in middle infrared * 

-3.30 

-3.60 

Rayleigh in near infrared * 

-3.12 

-2.82 


* Overlapping with water vapor absorption is included. 
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FIGURE CAPTIONS 

Figure 1 . Spectral distribution of the ozone transmission function of the entire atmospheric 
column in the vertical direction averaged over 10 cm' 1 . 

Figure 2. Same as Fig. 1 except for the water vapor transmission function. 

Figure 3. Same as Fig. 1 except for the 0 2 transmission function. 

Figure 4. Same as Fig. 1 except for the C0 2 transmission function. 
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